Shared genetic architecture between autoimmune disorders and B-cell acute lymphoblastic leukemia: insights from large-scale genome-wide cross-trait analysis

Background To study the shared genetic structure between autoimmune diseases and B-cell acute lymphoblastic leukemia (B-ALL) and identify the shared risk loci and genes and genetic mechanisms involved. Methods Based on large-scale genome-wide association study (GWAS) summary-level data sets, we observed genetic overlaps between autoimmune diseases and B-ALL, and cross-trait pleiotropic analysis was performed to detect shared pleiotropic loci and genes. A series of functional annotation and tissue-specific analysis were performed to determine the influence of pleiotropic genes. The heritability enrichment analysis was used to detect crucial immune cells and tissues. Finally, bidirectional Mendelian randomization (MR) methods were utilized to investigate the casual associations. Results Our research highlighted shared genetic mechanisms between seven autoimmune disorders and B-ALL. A total of 73 pleiotropic loci were identified at the genome-wide significance level (P < 5 × 10–8), 16 of which had strong evidence of colocalization. We demonstrated that several loci have been previously reported (e.g., 17q21) and discovered some novel loci (e.g., 10p12, 5p13). Further gene-level identified 194 unique pleiotropic genes, for example IKZF1, GATA3, IKZF3, GSDMB, and ORMDL3. Pathway analysis determined the key role of cellular response to cytokine stimulus, B cell activation, and JAK-STAT signaling pathways. SNP-level and gene-level tissue enrichment suggested that crucial role pleiotropic mechanisms involved in the spleen, whole blood, and EBV-transformed lymphocytes. Also, hyprcoloc and stratified LD score regression analyses revealed that B cells at different developmental stages may be involved in mechanisms shared between two different diseases. Finally, two-sample MR analysis determined causal effects of asthma and rheumatoid arthritis on B-ALL. Conclusions Our research proved shared genetic architecture between autoimmune disorders and B-ALL and shed light on the potential mechanism that might involve in. Supplementary Information The online version contains supplementary material available at 10.1186/s12916-024-03385-0.


Background
B-cell acute lymphoblastic leukemia (B-ALL) is a prevalent subtype of leukemia characterized by its highly malignant nature, primarily originating from the clonal expansion and abnormal proliferation of B lymphocytes within the hematopoietic system [1].Autoimmune disorders are characterized by a disruption in self-tolerance, resulting in pathological alterations and clinical symptoms arising from immune responses targeting self-components [2].Concurrently, the pathogenesis of several autoimmune disorders is intricately interwoven with the malfunctioning of B cells within the humoral immune system.The excessive activation of self-reactive B cells precipitates an overproduction of autoantibodies and immune complexes, which, in turn inflict damage upon a multitude of tissues and organs, culminating in the emergence of various autoimmune disorders [3].To summarize, B cells assume a pivotal role in the orchestration of humoral immune responses, and their deregulation markedly contributes to the onset of autoimmune diseases and B-cell malignancies [4].
Epidemiological investigations have discovered associations between autoimmune disorders and B-cell malignancies.For example, rheumatoid arthritis (RA) patients exhibit a twofold increased risk of concomitant B-cell lymphomas when compared to their healthy counterparts [5].In the case of systemic lupus erythematosus (SLE) and Sjögren's syndrome patients, the risk amplifies significantly to 2.7-7.5 times [6] and 9-18 times [6], respectively.Previous studies observed that the standardized incidence ratio of ALL was estimated to be 2.77 after RA onset [7].Studies also showed that at the time of diagnosis of malignancy, 15-30% of patients present with many of the typical features of rheumatic diseases [8].However, current research focused primarily on the onset of autoimmune diseases on hematological malignancies risk, particularly diffuse large B-cell lymphoma and follicular lymphoma.This leaves a clear gap in understanding the pleiotropic mechanisms and bidirectional causations between B-ALL (a disease also derived from B lymphocytes) and autoimmune diseases.Only Li et al. have reported the shared mechanism between autoimmunity and B-ALL, specifically demonstrating the essential role of DYRK1a in mediating the noncanonical NF-κB activation induced by BAFF [9].This underscores the existence of substantial knowledge gaps in this field, highlighting the urgent need to ascertain shared risk loci between these two disorders.It is worth noting that traditional clinical or epidemiological research may encounter challenges in ensuring the statistical effectiveness of such investigations.
Recently, the linkage disequilibrium (LD) score regression (LDSC) approach has been developed to indicate whether there exists a genetic correlation between the two types of disease [10].It is unclear whether the overall genetic correlation is attributable to a few loci or the entire genome.Few studies to date have systematically evaluated genetic overlap, shared susceptibility genes, and causality between autoimmune diseases and B-ALL.Cross-trait analyses that utilize the correlation of GWAS signals to study polyvalent genetic variants or loci between multiple traits have been shown to accurately identify shared loci between diseases or traits [11][12][13].These pleiotropic loci can be targeted for intervention to potentially prevent or treat these diseases simultaneously.Recently, a novel method called "PLACO" was developed to identify pleiotropy at the SNP-level based on a level-α intersection-union test (IUT) [14].Therefore, it is important to determine specific genetic variants or loci that lead to genome-wide genetic correlations or to delve into the shared genetic etiology of these two types of diseases.Our research flowchart is shown in Fig. 1.

Genetic overlap at the genome-wide level
We used LDSC to evaluate the genetic structure shared between autoimmune disorders and B-ALL [10].The LD scores used in LDSC were calculated based on genotypes of common SNPs from European ancestry samples in the 1000 genomes project [37].Standard errors (SE) were estimated by the jackknife method in LDSC which was further used to correct for attenuation bias.Intercept of LDSC was used to evaluate potential population overlaps between studies from different consortiums [10].It is worth noting that no actual population overlap between autoimmune disorders and B-ALL studies existed in our analysis.A likelihood-based method, called high-definition likelihood (HDL), can utilize GWAS summary statistics to estimate genetic associations, which could reduce the variance of genetic association estimates by about 60% compared with the LDSC method [38].
We further investigated whether SNP heritability of autoimmune diseases and B-ALL was enriched in specific cells and tissues using hierarchical LDSC regression.Stratified-LDSC (S-LDSC) was applied to different immune cell data to assess whether specific cell types had significant genetic enrichment in these tissues.We downloaded 54 human tissues datasets from GTEx [39] and 292 immune cell types from the ImmGen consortium [40] (including B cells, γ δ T cells, α β T cells, innate lymphocytes, myeloid cells, stromal cells, and stem cells).After adjusting for the baseline model and all gene sets, we assessed the significance of the SNP heritability enrichment estimated in each tissue and cell by using the regression coefficient Z-scores and corresponding P values.

Identification of pleiotropic loci and genes by using PLACO
A pleiotropic analysis under composite null hypothesis (PLACO) was used to identify pleiotropy among multiple autoimmune diseases and B-ALL at the SNP-level.SNPs reach genome-wide significant level (P < 5 × 10 -8 ) Fig. 1 Study workflow and were viewed as pleiotropic variants.The functional mapping and annotation (FUMA) of GWAS was used to determine the genomic regions of these risk variants (i.e., pleiotropic loci) [41].Also, a Bayesian colocalization analysis was conducted to determine the pleiotropic loci shared by autoimmune diseases and B-ALL [42].To explore the shared mechanisms of the identified loci, nearby genes were mapped based on lead SNPs within each locus.Also, a generalized geneset analysis of GWAS data (multi-marker analysis of genomic annotation, MAGMA) approach was used to determine the biological function of these pleiotropic loci.Specifically, we performed MAGMA gene analysis to identify pleiotropic genes by properly incorporating LD between markers and to detect multi-marker effects (P < 0.05/18,345 = 2.73 × 10 -6 ) [43].MAGMA gene-set analysis was performed to investigate the biofunction of lead SNPs [43], and a total of 10,678 gene sets including curated gene sets (c2.all) and go terms (c5.bp, c5.cc, and c5.mf ) from Molecular Signatures Database (MSigDB) were finally tested [44].Bonferroni correction was performed for all tested gene sets to avoid false positives (P < 0.05/10,678 = 4.68 × 10 -6 ).Metascape webtools (metascape.org)performed a pathway enrichment analysis to determine the function of mapped genes based on MSigDB [44].Genome-wide tissue-specific enrichment analysis was conducted based on 54 GTEx tissues [45] for the genome-wide pleiotropic results generated by PLACO.We also calculated the average expression (log 2 transformed) of all identified pleiotropic genes in all 54 GTEx tissues [45] and tested tissue specificity by differentially expressed genes (DEGs) in each tissue (up-and down-regulated DEGs were precomputed by the signs of the t-statistics).

Summary-based Mendelian randomization
Summary-based Mendelian randomization (SMR) [46] method combined summary-level data from GWAS with data from expression quantitative trait loci (eQTL) studies to identify genes whose expression levels are associated with complex traits due to pleiotropy.It employs SMR and HEIDI methods to test pleiotropic associations between gene expression levels and complex traits of interest using summary-level data from GWAS and eQTL studies.This approach could be interpreted as an analysis to test whether the magnitude of SNP effects on phenotype is mediated by gene expression.

Multi-trait colocalization analysis
We utilized hypothesis prioritization for multi-trait colocalization (HyPrColoc) [47] method to perform multi-trait colocalization analysis to pinpoint the crucial roles that immune traits played in the onset of autoimmune disorders and B-ALL.Immune-wide GWAS data contains a total of 731 immune cells [48], which could be publicly available from the GWAS catalog (GCST0001391 ~ GCST0002121).Detailed information on the GWAS summary datasets for immune cells was added to Additional file 1: Supplementary Methods.

Causal association analysis
We performed a one-directional two-sample Mendelian randomization (MR) analysis to assess possible causal effects of autoimmune disorders on B-ALL risk.The "clumping" procedure in PLINK 1.9 software was used to extract independent significance SNPs for all autoimmune diseases (P < 5 × 10 -8 ), where r 2 was set to 0.001 and window size was set to a physical distance of 10,000 KB [49].Notably, r 2 was calculated based on the 1000 genomes project phase 3 as a reference panel.Proportion of variance explained (PVE) and F statistic (F > 10) was used to measure the strength of instrumental variables (IVs) (see Additional File 1: Supplementary Methods) [50].To verify causality among these trait pairs, six MR approaches were performed with each set of IVs, i.e., IVW, Debiased-IVW (DIVW) [51], weighted median approach [52], MR pleiotropy residual sum and outlier (MR-PRESSO) [53], MR-Egger [54], MR robust adjusted profile score (MR-RAPS) [55], and mode-based estimate [56] method.Cochran's Q statistics was used to examine the effect size heterogeneity across the IVs (see Additional File 1: Supplementary Methods) [57,58].Additionally, the intercept of MR-Egger regression and global test of MR-PRSSO were utilized to detect horizontal pleiotropy [53,54].Detailed information on used MR methods was described in Additional file 1: Supplementary Method.

Shared genetic architecture between autoimmune disorders and B-ALL
We first evaluated the genetic correlation between autoimmune diseases and B-ALL and results from the LDSC and HDL methods were highly consistent (Table 1 and Additional file 2: Table S2).Specifically, by using the LDSC method, six traits were identified to be genetically correlated with B-ALL, including AOA, HT, IBD, CD, RA, and MS.While implementing the HDL method, significant genetic correlations were observed among AOA, HT, PBC, RA, MS, and B-ALL, leading to a final union set of seven pairwise traits for further analysis.However, we did not find significant genetic correlation between IBD and CD and HDL results.It was noting that only RA remained significantly genetical correlated with B-ALL risk after applying the Bonferroni correction (P = 0.003 < 0.05/16).

Pleiotropic loci and genes identified for multiple autoimmune disorders and B-ALL
Given the shared genetic mechanisms between autoimmune diseases and B-ALL identified by LDSC and HDL, we used novel pleiotropy analyses (PLACO) to identify potential pleiotropic loci for both diseases (Additional file 1: Fig. S1).The QQ plots demonstrated no premature divergence between observed and expected values, ruling out the possibility of group stratification (Additional file 1: Fig. S2).Based on PLACO results, we identified a total of 73 pleiotropic genomic risk loci associated with both autoimmune disorders and B-ALL using FUMA (P < 5 × 10 -8 ) (Fig. 2, Additional file 1: Fig. S1, and Additional file 2: Table S3).Colocalization analysis finally identified 16 of 73 (21.9%) potential pleiotropic loci with PP.H4 greater than 0.7 (e.g., 5p13) (Table 2).The regional plots for each trait pair are presented in Additional file 1: Fig. S3 ~ S8.Notably, some pleiotropic regions were shared between different pairs, for example, genome regions 7p12.2, 10p14, 6q27, and 10p12.31were identified in four pairs (Additional file 2: Table S4).The MAGMA analysis of gene sets suggested that the identified pleiotropic loci may participate in the control of the immune system, hematopoiesis, and various other processes (Fig. 3A and Additional file 2: Table S5).Notably, significant monocyte differentiation pathway was found for all trait pairs, and significant leukocyte differentiation was found for all five trait pairs.Further tissuespecific analysis found these risk loci were enriched in several immune-related tissues (e.g., spleen, whole blood, Epstein-Barr virus (EBV)-transformed lymphocytes) (Fig. 3B and Additional file 2: Table S6).ANNOVAR category annotation showed that 28 of 73 lead SNPs (38.4%) were intronic variants and 30 of 73 (41.1%) were intergenic variants.Only 2 of 73 (3%) lead SNPs were exonic variants (Additional file 2: Table S3).

Pleiotropic genes identified and enrichment analysis
We used different methods to map the identified SNP-level signals into the gene-level signals.By using MAGMA gene analysis, a total of 341 significant pleiotropic genes were determined as pleiotropic genes between multiple autoimmune diseases and B-ALL (194 unique) (Additional file 2: Table S7 and Additional file 1: Fig. S9).Additional file 2: Table S8 lists the details of these genes.MAGMA gene analysis detected 92 repeated pleiotropic genes across different trait pairs, with IKZF1 identified as a pleiotropic gene for six pairs, followed by MLLT10, FIGNL1, RNASET2, CCR6, GATA3, CLN3, PIP4K2A, DDC, RP11-514O12.4,FGFR1OP, and GRB10 in four trait pairs.eQTL analysis identified multiple hits of IKZF1 in blood-and immune-related tissues (e.g., naïve B cell, CD19 B-cell, EBV-transformed lymphocytes cells, cis-eQTLs, trans-eQTLs, spleen, whole blood).Five genes (i.e., TUFM, ZC2HC1A, RNASET2, GSDMB, and ORMDL3) were observed to be significant in five different tissues.We summarized the landscape of pleiotropic genes identified in different methods and tissues in Fig. 4. We observed several genes (RNASET2 and FIGNL1) were significantly mapped in different tissues with different methods.The IKZF1 gene was also highlighted in whole blood tissues (Fig. 4).
The shared mechanism between autoimmune diseases and B-ALL may involve specific organs or tissues involvement.Numerous genes (e.g., TOP2A, IKZF3, MYB, and CD80) showed significant differential expression in EBV-transformed lymphocytes, and APOBR, IKZF1, and IL7R showed significant differential expression in spleen and whole blood tissues (Additional file 1: Fig. S10 and Additional file 2: Table S9).Tissue enrichment analysis showed that these genes were also enriched into the spleen and EBV-transformed lymphocytes (Additional file 1: Fig. S11 and Additional file 2: Table S10).Additional S-LDSC based on multiple tissues identified significant SNP heritability enrichment for all autoimmune diseases (except AOA) in each of the monocytes, blood cells, and spleen, after adjusting for the baseline Fig. 2 The circular diagram presents pleiotropic loci and genes identified by PLACO among seven trait pairs.Note: Shared loci identified by colocalization analysis are highlighted in orange; shared genes identified by MAGMA analysis are highlighted in blue.B-ALL B-cell acute lymphoblastic leukemia, AOA adult-onset asthma, HT hypothyroidism, PBC primary biliary cirrhosis, IBD inflammatory bowel disease, CD Crohn's disease, RA rheumatoid arthritis, MS multiple sclerosis model (Additional file 1: Fig. S12 and Additional file 2: Table S11).Further enrichment analysis of the GO biological processes associated with these genes indicated higher enrichment in the cellular response to cytokine stimulation, B cell activation, response to tumor necrosis factor, inflammatory response, and receptor signaling Table 2 16 colocalized loci identified by colocalization analysis performed on 73 pleiotropic loci (PP.H4 > 0.7) Lead SNP was the SNP with minimum P values within the corresponding locus.PP.H4 was the posterior probability of H4 calculated by coloc analysis; the Locus boundary was defined as "chromosome: start-end" PP.H4 the posterior probability of H4, B-ALL B-cell acute lymphoblastic leukemia, AOA adult-onset asthma, HT hypothyroidism, PBC primary biliary cirrhosis, IBD inflammatory bowel disease, CD Crohn's disease, RA rheumatoid arthritis, MS multiple sclerosis  pathway via JAK-STAT (Fig. 5A).These pathways play important roles in immune regulation and leukemogenesis.Cell type enrichment analysis showed the highest significance for bone marrow naïve T cells (Fig. 5B).Furthermore, we found that these genes were numerically enriched in several immunologic signatures (e.g., MEM-ORY VS CD21HIGH TRANSITIONAL BCELL DN) (Fig. 5C).The PPI analysis showed that five PPI networks were constructed, including the JAK-STAT signaling pathway and multiple pathways related to DNA damage were involved.And 22 proteins (e.g., STAT, NFKB1, and GATA3) could participate in these pathways (Fig. 5D).Also, the results suggest that heritability is enriched in the blood, EBV-transformed lymphocytes, whole blood, and palatine tonsil tissues among five or more autoimmune diseases and B-ALL.

Immune-related mechanisms shared between autoimmune disorders and B-ALL
The shared mechanism involving affected tissues such as the spleen, lymphocytes, and whole blood, suggested an important involvement of immune mechanisms in the inter-disease.We used the S-LDSC method to determine the heritability enrichment of pleiotropy in immune cells and the HyPrColoc method to identify immune cells with co-localization signals with pleiotropic motifs.S-LDSC observed heritability enrichment of B cells in both autoimmune diseases and B-ALL.When analyzing the enrichment of immune traits from ImmGen, we also observed that two cell traits in the B cell panel were enriched: B.FrE.BM (CD19 + IgM + AA4.1 + HSA + ) and preB.FrD.BM (CD19 + IgM − CD45R + CD43 − ).Additionally, numerous cell traits in the T cell panel were also identified, implying the potential immune mechanisms shared (Additional file 1: Fig. S12 and Additional file 2: Table S11).Then multi-trait colocalization analysis by using HyPrColoc was performed to pinpoint key immune cells (Additional file 2: Table S12).

The causal relationship between autoimmune diseases and B-ALL estimated by MR
MR analyses using the IVW method showed significant positive associations between two autoimmune diseases (AOA and RA) and B-ALL risk (Fig. 6A and Additional file 2: Table S13).The risk of B-ALL was found to be able to be increased as the risk of AOA increases, the effect size was estimated by using the IVW method (OR = 1.223, 95%CI = 1.048 ~ 1.426, P = 0.010).Another four methods (DIVW, MR-RAPS, MR-PRESSO, and slope of MR-Egger) are consistent with the results of the IVW method.Although a significant intercept of MR-Egger might indicate the existence of potential horizontal pleiotropy, the global test of MR-PRESSO ruled out this possibility (P = 0.632).We also observed significant causal effects of RA onset on B-ALL risk by using  S14).Additional scatter and funnel plots eliminate the possibility of potential outliers (Fig. 6C-D).However, after the Bonferroni adjustment, no causal associations between autoimmune disorders and B-ALL remained statistically significant (P = 0.003 < 0.05/16).Finally, reverse MR analysis ruled out the possibility of reversedirectional causality.

Discussion
Given the critical contribution of B cells to autoimmune disorders and B-ALL, there may be a complex relationship between them [60].The study employed comprehensive genetic methodologies to investigate the genetic correlation between autoimmune disorders and B-ALL.
The study determined pleiotropic loci using cross-trait PLACO analysis and identified pleiotropic genes through the MAGMA method.Then the key pathways and immunological mechanisms involved were identified.Finally, comprehensive MR analysis and sensitive analysis established the causal relationships between autoimmune diseases and B-ALL.
Through genetic correlation analysis, we observed significant genetic overlap between B-ALL and seven autoimmune disorders, including AOA, HT, IBD, CD, PBC, RA, and MS.We provide strong evidence for a shared genetic mechanism between RA and B-ALL, as well as MR evidence suggesting that patients with RA symptoms should be alerted to the risk of progression to ALL, which is consistent with previous studies [7,8].Additionally, study have shown that 34 of 699 ALL patients diagnosed and followed had previously received varying doses of steroids for aplastic events or arthritis-based rheumatic diseases [61].By using genetic variables, MR methods could well avoid the influence of possible confounding factors.Therefore, we believe that in addition to the effect of immunosuppressants, RA itself will also play an important role in the risk of B-ALL.We also observed significantly causal effects of AOA on B-ALL risk, which was ambiguous in previous studies: a systematic review supported the protective effect of asthma on ALL [62], two types of research included showed significant high risks of ALL in patients with a history of asthma [63,64].
We identified a series of genetically risk loci associated with both autoimmune diseases and B-ALL, and some of which were observed in multiple phenotype pairs (e.g., 7p12.2, 10p14, 6q27, 10p12.31).Previous studies gave the evidence of key role these loci played in the development Fig. 6 A The forest plot shows causal associations between autoimmune disorders and B-ALL by using one-directional MR analysis.Note: Causal effects were estimated by using IVW method.B Forest plot shows causal effects of AOA and RA on B-ALL risk estimated by using different methods.C Scatter plot shows significant causal association between AOA and B-ALL risk.D Funnel plot shows significant causal association between AOA and B-ALL risk.E Scatter plot shows significant causal association between RA and B-ALL risk.F Funnel plot shows significant causal association between RA and B-ALL risk.Associations highlighted with red represent that associations were significant in more than three main MR methods of autoimmune disorders and B-ALL.For example, loci on 7p12.2 (IKZF1) had been proven to be associated with risk of childhood B-ALL [65], which was also identified as susceptibility genes for SLE [66].After searching for the GWAS catalog, 7p12.2 had been reported to be associated with multiple autoimmune disorders, including CD [23], IBD [23], RA [25], and MS [67].GATA3 (10p14) is a key regulator in the immune system, especially in the differentiation and function of type 2 helper (Th2) cells [68].Th2 cells have been demonstrated to play a role in various autoimmune diseases, including SLE and IBD [69,70].Recent research also highlighted the role of noncoding genetic variation (rs3824462) in GATA3, linking it to an increased risk of Ph-like ALL, a common subtype of B-ALL.The study revealed that rs3824462 induced local and global changes in chromatin conformation, activating JAK-STAT pathway and promoting disease development [71].
We searched for the identified risk loci in the GWAS catalog (last update in 2023 December 20) [72] and found that some of the risk loci have been reported to be associated with both B-ALL and autoimmune disorders (Additional file 1: Fig. S13 and Additional file 2: Table S15).For instance, the 17q21 locus is implicated in various autoimmune diseases, including asthma [73,74], IBD [75,76], T1D [77], and SLE [78].This locus, housing IKZF3, GSDMB, and ORMDL3, involved in lymphocyte development [79], pyroptosis [80], and inflammatory response [81], has been challenging to dissect.GSDMB and ORMDL3 represent the target genes of rs2290400, and its minor allele is associated with a protective effect against ALL [82].IKZF3 polymorphism contributes to B-ALL with a 1.5-fold to twofold increase in relative risk [83].Genes previously reported to be associated with leukemia have also been observed in our results to be correlated with autoimmune diseases: MLLT10 (10p12) participates in various chromosomal rearrangements associated with ALL and acute myeloid leukemia (AML) [84].It is implicated in chromatin structure regulation and DNA damage response, deemed crucial for early development, maintenance, and differentiation of hematopoietic stem cells.While direct evidence for the impact of MLLT10 on autoimmune diseases has not been established, studies indicated a close association with C-reactive protein levels [85], widely recognized as a valuable indicator of disease activity in various autoimmune rheumatic diseases [86].Simultaneously, certain genes previously reported to be associated with autoimmune disorders have also been found in our results to be associated with B-ALL.IRGM (5q33) encodes a member of the immunity-related GTPase family, crucial in innate immunity and inflammatory responses [87].Previous studies have linked IRGM to CD [88][89][90], UC [91], and IBD [23].CAPSL (5p13) has been reported to be associated with PBC [92], T1D [93], asthma [94], and SLE [95].Although direct evidence of its association with ALL is lacking, increased mRNA levels have been observed in AML patients [96].Additionally, the long non-coding RNA C5orf56 (5q31) has been identified for its protective role in IBD [97].SCHIP1 (3q25) has been associated with SLE [98], while RNASET2(6q27) has been identified as a risk gene for both vitiligo [99] and GD [100].
Shared genetic structures observed in our research revealed common mechanisms between autoimmune disease and B-ALL.Identified genes were observed to participate in several pathways, like B cell activation, cellular response to cytokine stimulus, and inflammatory response.For each disease pair, we observed a significant enrichment of pleiotropy to the spleen, a critical site for B cell development.Notably, a substantial presence of BAFF-R-associated immune signature, a key regulator of B cell function and survival, was discerned in a multitrait colocalization analysis.These findings collectively underscored the pivotal role played by B cells in both autoimmune disorders and B-ALL.In autoimmune conditions, B cells are exposed to antigens, undergo activation, and subsequently proliferate and expand clonally, thereby increasing the risk of accumulating genetic mutations, and finally leads to the emergence and progression of B-ALL [60].We can think that ORMDL3 and IKZF3, mentioned earlier, play crucial roles in this context, as evidenced by prior literature reporting ORMDL3's vital role in B cell survival [101], and IKZF3's predominant regulation of B cell differentiation, activation response, and proliferation [102].Furthermore, malignancies arising from B cells consistently exhibit concurrent autoimmune disorders at any stage, whereas those derived from T cells are less commonly linked to autoimmune phenomena [103].Nevertheless, our findings also identified numerous cell traits in the T cell panel, and we speculate that this may be attributed to interactions between B and T cells.The JAK-STAT pathway may represent a crucial mechanism in this context, as it has been targeted in autoimmune diseases [104] and its role in B-ALL involves the disruption of preleukemic cells differentiation [105].Our results highlighted the critical role of EBV infection as a trigger for both autoimmune disorders and B-ALL: tissue-specific analysis revealed enriched risk loci in EBV-transformed lymphocytes, and the central role of IKZF1 in this cell was also identified by gene-level analyses.EBV remains latent in memory B cells after infection, and reactivation can induce B cell clonal immortalization, promoting lymphomagenesis [106].Additionally, EBVinduced autoimmunity has been reported to increase the risk of autoimmune diseases [107].

Limitations
Our study is not without limitations.Firstly, as with other similar studies, the data used in this study was summarylevel, and individual-level datasets were not available.Further stratification of the population (e.g., gender, age, etc.) was therefore not possible.Secondly, the sample size of immune cell GWAS used in this study was limited.Therefore, caution should be exercised in interpreting the role of immune cells and drawing conclusions in our studies.Thirdly, it should be noted that our study was limited to European ancestry and may not be generalizable to other ancestries.It is important to be equally cautious in concluding our findings since the relatively small sample size of B-ALL may result in limited statistical power.

Conclusions
Our research has uncovered the intricate connections between autoimmune disorders, especially AOA, HT, IBD, CD, RA, and MS and B-ALL.Identification of pleiotropic risk loci (7p12, 10p14, 6q27, and 10p12) and genes (IKZF1, GATA3, IKZF3, GSDMB, and ORMDL3) shared between diseases suggested shared mechanisms, such as B cell activation and JAK-STAT pathway, common triggers like EBV infection.Additionally, our findings have shed light on and the causal links between autoimmune disorders (AOA and RA) and B-ALL.

Fig. 3
Fig. 3 Bar plot of MAGMA gene-set (A) and tissue-specific (B) analysis for genome-wide pleiotropic results.Note: The red dotted line represents the significance of 0.05 after multiple corrections, and the blue represents the significance of 0.05.B-ALL B-cell acute lymphoblastic leukemia, AOA adult-onset asthma, HT hypothyroidism, PBC primary biliary cirrhosis, IBD inflammatory bowel disease, CD Crohn's disease, RA rheumatoid arthritis, MS multiple sclerosis

Fig. 4
Fig. 4 Overview of pleiotropic genes (highlighted in all three signals) for the autoimmune disorders and B-ALL.Note: The signals represent hits of genes across different trait pairs.eQTL expression quantitative trait loci, SMR summary-based Mendelian randomization, AD autoimmune disorders, B-ALL B-cell acute lymphoblastic leukemia, AOA adult-onset asthma, HT hypothyroidism, PBC primary biliary cirrhosis, IBD inflammatory bowel disease, CD Crohn's disease, RA rheumatoid arthritis, MS multiple sclerosis Results highlight 59 pleiotropic loci, of which 19 were unique, and these loci support the important role of 42 unique immune cells in autoimmune diseases and B-ALL by sharing causal variants.Our results support the critical influence of BAFF-R, CD4, CD45, and CD28 on different cells.Notably, a total of six BAFF-R-related immune traits were observed, including BAFF-R on B cell, BAFF-R on CD20 − , BAFF-R on CD24 + CD27 + , BAFF-R on IgD + CD24 − , BAFF-R on IgD + CD24 + , and BAFF-R on IgD + CD38 − .Interestingly, BAFF-R on B cell and BAFF-R on CD24 + CD27 + were both shared among three trait pairs (i.e., B-ALL&IBD, B-ALL&PBC, B-ALL&RA).

Fig. 5 A
Fig.5A Pathway enrichments for identified pleiotropic genes (KEGG, GO, Wiki pathways).B Cell-type enrichments for identified pleiotropic genes.C Immune signatures enrichments for identified pleiotropic genes.D Protein-protein interaction analysis based on identified pleiotropic genes